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We investigate the stability of self-similar solutions for a gravitationally col- 
lapsing isothermal sphere in Newtonian gravity by means of a normal mode 
analysis. It is found that the Hunter series of solutions are highly unstable, 
while neither the Larson-Penston solution nor the homogeneous collapse one 
have an analytic unstable mode. Since the homogeneous collapse solution is 
known to suffer the kink instability, the present result and recent numerical 
simulations strongly support a proposition that the Larson-Penston solution 
will be realized in astrophysical situations. It is also found that the Hunter 
(A) solution has a single unstable mode, which implies that it is a critical so- 
lution associated with some critical phenomena which are analogous to those 
in general relativity. The critical exponent 7 is calculated as 7 ~ 0.10567. In 
contrast to the general relativistic case, the order parameter will be the col- 
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lapsed mass. In order to obtain a complete picture of the Newtonian critical 

phenomena, full numerical simulations will be needed. 
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I. INTRODUCTION 



Spherically symmetric self-similar systems have been widely studied in the context of 
both Newtonian gravity and general relativity. Self-similar solutions in Newtonian gravity 
have been studied in an effort to obtain realistic solutions of gravitational collapse leading to 
star formation In this context, Larson and Penston independently found a self-similar 

solution, which describes a gravitationally collapsing isothermal gas sphere Thereafter, 
Hunter found a new series of self-similar solutions, and that a set of such solutions is infinite 
and discrete |Q. Whitworth and Summers investigated the mathematical structure of the 
equation for self-similar solutions in more detail and found a band structure of a new family 
of self-similar solutions with loss of analyticity |^. The solutions were classified into two 
types, based on the behavior around the sonic point. These self-similar solutions were 
generalized to general relativity by Ori and Piran 

In general relativity, self-similar solutions called attention in the discovery of the critical 
behavior by Choptuik P] . Evans and Coleman found similar critical behavior in the collapse 
of a radiation fluid 0. A renormalization group approach showed that the critical solution, 
which is at the threshold of the collapse to a black hole, is an analytic self-similar solution 
with the unique unstable mode and that the critical exponent which appears in the scaling 
law of the formed black hole mass is equal to the inverse of the eigenvalue of the unstable 
mode for a perfect fluid case ||10|-Q . In the previous work |13| , the authors found that the 
general relativistic counterpart of the Hunter (A) solution has a single unstable mode and 
expected that the Hunter (A) solution is a critical solution of the Newtonian counterpart of 
the critical behavior in gravitational collapse. The first purpose of this paper is to confirm 
this expectation and to calculate the accurate value of the critical exponent. 

In addition to this critical nature, we should mention the role of a self-similar solution 
as an attractor of gravitational collapse. Recent numerical simulations and results of mode 
analyses showed that the Larson-Penston solution is the best description for the central part 
of a collapsing gas sphere in Newtonian gravity [p^ -p!8[] and in general relativity |13|,|19 
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Hanawa and Nakayama examined the spherically symmetric unstable modes in linear 



order of which the growth is faster than as t ^ 0. They showed that the Larson- 
Penston solution has no such mode, while the Hunter (B) and (D) solutions have some 
unstable modes by means of a normal mode analysis. They concluded that the Hunter (B) 
and (D) solutions are unstable and not likely to be realized in generic situation, while the 
Hunter (A) and (C) solutions were dropped from the analysis. There is no reason to rule out 
the Hunter (A) and (C) solutions a priori . The second purpose of this paper is to complete 
the stability analysis on the self-similar solutions including these solutions. 

The organization of this paper is the following. In section II, basic equations are pre- 
sented. In section HI, eigenvalues of unstable modes for the self-similar solutions are pre- 
sented. Section IV is devoted to discussions. In Section V, we summarize the paper. 

II. BASIC EQUATIONS 
A. Basic equations in the zooming coordinates 

A gravitationally collapsing isothermal sphere is described in spherical coordinates by 

| + ^|(^»-0. (2.1) 

|(p«) + 4#('-V«Vc2|^ + P^ = 0. (2.2) 

ot or or 

dM DM , , 

^ = fcrV, (2.4) 

where p, v, M, Cg and G denote the density, radial velocity, total mass inside the radial coor- 
dinate r, sound speed, and gravitational constant, respectively. We introduce the zooming 
coordinate z = {cst)/r. We also introduce dimensionless functions U,P and m: 

vir,t) = -cMir,t), (2.5) 
P(M)^43^. (2.6) 
M(nt) = *^. (2.7) 
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Then, the equations are expressed as 

+ {1 + zU)P' + zPU' = 0, (2.8) 
-t{UP + UP) - z{U + zU"^ + z)P' 

-zP{l + 2zU)U' -2zP + mPz^ = 0, (2.9) 

tm + zm' + m- PU = 0, (2.10) 

-z^m' = P, (2.11) 

where the prime and dot denote the partial derivatives with respect to z and t, respectively. 

B. Self-similar solutions 

Self-similar solutions are characterized hj U = U{z),P = P{z) and m = mlz). From 
this ansatz, equations (p.8|)-( pn| ) become 



jzU + l)[P{zU + 1) - 2] 

^ ~ (zu + iy-z^ ' ^ ^ 

zP[2 - P{zU + 1)] 
{zU + 1)2 - z^ 

m = P{U +l/z). (2.14) 



f'- r (2.13) 



If the analyticity at the center is assumed, self-similar solutions can be expanded in Taylor 
series around z —>■ —oo, i.e., t — oo or r ^ as follows: 

3z 45 V3 J z^ \z^J' 
Q = Hz^P) = Qo + l{l- e«") ^ + O (^) . (2.15) 

Therefore, self-similar solutions which have regular center are specified by one parameter 
Qq. Taylor series expandability in the neighborhood of the sonic point z = Zg {zU + z = — 1) 
requires that there are two possible analytic solutions. The first solution (type 1) is 



and the second one (type 2) is 

U = - + - + (^s - Zs - l){z - 

\ Zs) Zs 2^3 (3z, + 2) 

p 2 2{z-z,) (7zl + Qzs + l){z-z,f 

Zs zl (3^, + 2) • ^ • ^ 

Therefore, self-similar solutions which are analytic at the sonic point are specified by one 
parameter Zs around the sonic point. We can find one exact solution and numerical ones 
with analyticity both at the center and at the sonic point. The former is a homogeneous 
collapse solution: 

Qo = ln^, z, = -\ P=^. f^ = -|-' "^ = A- (2-18) 

3 3 6z^ 3z 9z-^ 

The values of Qo and Zg for numerical solutions (the Larson- Penston solution and the Hunter 
(A)-(D) solutions) are summarized in Table |T[ These self-similar solutions are displayed in 
Figs. 

The homogeneous collapse solutions is the only solution which has the big crunch sin- 
gularity. The big crunch occurs at t = 0, i.e., the singularity occurs at the same time 
everywhere. Unlike the homogeneous solution, the Larson-Penston solution and the Hunter 
(A)-(D) solutions are regular at t = 0, except for at r = 0. The Hunter (A) and (C) solutions 
encounter another sonic point. Neither the Hunter (A) nor (C) solutions can pass through 
the second sonic point regularly (see Figs. |l] and §). The Hunter (A)-(D) solutions are char- 
acterized by the number of oscillations in their profiles zU (see Fig. 0). Unlike the Hunter 
series, no oscillation can be seen in the profile for the Larson-Penston solution. This means 
that the Larson-Penston solution is a pure collapse solution. A similar property is shown 
also in P (see Fig. ^). In the Larson-Penston solution and the Hunter (A)-(D) solutions, 
the density profile has the maximum value at the center and decreases monotonically with 
increasing in z (see Fig. ^). In these solutions, the Hunter (D) solution has the biggest 
central value of P. We note that the Hunter solutions and the homogeneous collapse one 
are type 1, while the Larson-Penston solution is type 2. 
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C. Perturbation equations 



We consider the spherically symmetric linear perturbations around the self-similar solu- 
tion. We define the perturbation quantities as 

P{r,t) = Po{z)+eP,{t,z) + 0{e^), 
U{r, t) = Uo{z) + eUi{t, z) + 0{e''), 

m{r, t) = mo{z) + emi{t, z) + 0{e^), (2.19) 

where Pq, Uq and mo are the background self-similar solution and £ is a small parameter 
which controls the expansion. Then we find the equations for perturbations up to linear 
order of £ as 

tP 

— + zUiP' + P[{1 + zU) + zPU[ + zPiU' - 0, (2.20) 
z 

-zP'(Ui + 2zUUi) - zP[{U + zU'^ + z) 
-tUiP - tUPi - z{PiU' + U[P){1 + 2zU) 

-2z^UiPU' - 2zPi + z^(Pmi + mPi) = 0, (2.21) 

zm[ + mi + tnii - {UPi + PUi) = 0, (2.22) 

-z'^m[ = Pi, (2.23) 

where we have omitted the suffix o for simplicity. 

We assume the time dependence of the perturbations as 

Pi{r,t) = SP{z)e''\ 
Ui{r,t) = SU{z)e''\ 

mi{r,t) = 5m{z)e''^, (2.24) 
where r = — ln(— t). Then we find the following equations for the perturbations: 



[(1 + zUf - z^]SP' 



2z - Pz{l + zU) + -(1 + zU) - 7^(1 + zU) 

Z 1 — (7 



5P 



+ 



-aP + z^PU' - (1 + zU)zP' - 



fP^ 
1 - a 



6U, 



(2.25) 



-zP[{l + zU)"^ - z'^]6U' 



[1 + zUf (U'-P 



1 - a 



2(1 + zU)+a- z^U' 



6P 



+ 



[1 + zU) -aP + z^PU' - 



IP^ 
I- a 



z'P' 



6U, 



-z^5m' = 6P, 

(1 - a)6m = Q + f/^ 



6P + P6U. 



(2.26) 
(2.27) 
(2.28) 



Here we examine boundary conditions which the perturbations should satisfy at the 
boundaries. First we consider the regular center [z = — oo). The perturbations near the 
center must satisfy 



SU 
6P 



az^ 



(2.29) 



Next we consider the sonic point {z = Zs). We require that the perturbation of the density 
gradient is finite. The boundary condition for the perturbations at the sonic point is 



z 



6P + 



{l-a){-z^P' + aP-z^PU') 



+ P' 



su = o. 



(2.30) 



Only for a discrete set of a, there exists a solution of perturbation equations which is regular 
both at the regular center and at the sonic point. Thus we can obtain eigenvalues a and 
the associated eigenmodes. It is easily shown that the homogeneous collapse solution has 
one stable mode a = —2/3 (see Appendix A). We note that all self-similar solutions have a 
ghost mode with a = 1 (see Appendix B). 



III. NUMERICAL RESULTS 



We have solved equations ( |2.25| ) and ( |2.26| ) numerically by fourth-order Runge-Kutta-Gill 
integration with adaptive stepsize control and obtained eigenvalues of them by the shooting 
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method. We have set SUq = 1. When the amphtude of the perturbations has grown up to 
larger than lO^*', we have appropriately scaled down their amplitude equally in order for the 
calculation not to be stopped by overflowing error. We have used equation ( p.28| ) to check 
the numerical error. The numerical error of this calculation has been within 10~^. We have 



assumed that the eigenvalue is positive. Since a has upper bound Wexp((5o) + 1 for the self 



similar solutions ||T5|, we have sought the eigenvalues in the region < cr < ^Jexp{Qo) + 1. 
The result is summarized in Table |T|. 
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TABLES 


Solution 


Qo 


v/oxp(Qo) + 1 


'-S 


Homogeneous 


In (2/3) 


1.8165 


-1/3 


Larson-Penston 


0.5101 


2.2905 


-0.4271 


Hunter (A) 


7.45616 


42.599 


-1.35305 


Hunter (B) 


11.236 


276.34 


-0.9071 


Hunter (C) 


16.322 


3502.7 


-1.0295 


Hunter (D) 


20.975 


35865 


-0.9911 




TABLE L Qo and 


Zs of the self-similar solutions. 




Solution 




Mode 


a 


Homogeneous 






Nothing 


Larson-Penston 






Nothing 


Hunter (A) 




1 


9.4637 


Hunter (B) 




1 


5.49 






2 


5.74 X 10^ 


Hunter (C) 




1 


6.56 






2 


5.89 X 10^ 






3 


7.18 X 10^ 


Hunter (D) 




1 


6.22 






2 


5.85 X 10^ 






3 


5.95 X 10^ 






4 


7.35 X 10^ 



TABLE n. Summary of stability analysis (unstable modes). 
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It is found that neither the homogeneous collapse solution nor the Larson-Penston one 
have no unstable mode and that each Hunter solution has the same number of unstable 
modes as the number of oscillations in their profiles of zU. We obtain the eigenvalues up 
to five digits only for the Hunter (A) solution since this relates to the critical exponent in 
critical phenomena. We number the unstable modes for each Hunter solution in order of 
magnitude of their eigenvalue. The eigenvalues of mode n{n = 1,2,3) for the Hunter (B), 
(C) and (D) solutions are, if there are, close to each other. It can be seen that the eigenvalue 
of mode m(m = 2,3,4) is, if there is, approximately ten times greater than the eigenvalue 
of mode (m — 1) for the Hunter (B), (C) and (D) solutions. 

The mode functions of the perturbations for the Hunter (A) solution are shown in Figs. 
H and 1^. The mode function of the density perturbation has a node. It has a large amplitude 
near the center and a small amplitude with the opposite sign near the sonic point. This 
perturbation makes the concentration of the density strong or weak. The mode function of 
the velocity perturbation does not have a node. It is found that the "positive" perturbation 
enhances the collapse of the gas spheres, while the "negative" one promotes the gas to 
disperse away. 



IV. DISCUSSIONS 

We have investigated the stability of the Larson-Penston solution, the homogeneous 
collapse solution and the Hunter (A)-(D) solutions by a normal mode analysis, assuming 
that the eigenvalue is positive. This assumption can be justified for Sfto" > 1 (See Appendix 
in [|15|). In addition, the results of the Lyapunov analysis by [jl2[ for Cg <^ c, where c is the 



speed of light, strongly suggests the validity of this assumption at least for the Hunter (A) 
solution. 

From the results in the previous section, the Hunter (A)-(D) solutions are unstable and 
not likely to be realized. It has been shown that both the homogeneous collapse solution 
and the Larson-Penston solution have no unstable mode. Since the homogeneous collapse 
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solution has kink instability as proved by Ori and Piran pO[ , the Larson- Penston solution is 
the best description of the central part of generic spherical collapse of isothermal gas. This 
result is consistent with the numerical simulations P JI^ , p!6[] . 

The eigenvalues of the Hunter (B) and (D) solutions obtained here are slightly different 
from the result of Hanawa and Nakayama |T^ , however the eigenvalues of the Hunter (A) and 



(B) solutions obtained here are consistent with those of the Newtonian limit in the general 
relativistic analysis |[Tl|JT3|JT^ . We believe that the eigenvalues obtained in the present paper 



are quite accurate. 

Whitworth and Summers claimed that the Hunter (A) solution (and also the Hunter (C) 
solution) is unacceptable since it cannot pass through the second sonic point regularly (see 
Figs. H and 0) . However, we note that we can prepare regular initial density profile which 
develops these solutions. The fact that these solutions cannot pass thorough the second 
sonic point analytically only implies that the self-similarity and regularity requirements are 
incompatible in the evolution of the sphere for t > and does not rule out the possibility that 
they may describe generic collapse of the isothermal sphere leading to the core formation. 
Only by the results of the present analysis, we can deny this possibility. 

We have found that the Hunter (A) solution has a single unstable mode. From the 



discussions in use of a renormalization group [|10l, it implies that this solution is a critical 



solution of some critical phenomena. The obtained eigenvalue of the unstable mode for the 
Hunter (A) solution is 9.4637, from which the critical exponent is calculated as 7 = a^^ ~ 
0.10567. The order parameter in this Newtonian case will be the collapsed mass in place of 
the formed black hole mass in general relativity. It is clear that full numerical simulations 
will give a complete picture of the Newtonian critical behavior. 

V. SUMMARY 

It is shown by means of a normal mode analysis that the Hunter (A)-(D) solutions are 
unstable, while neither the Larson-Penston solution nor the homogeneous collapse one have 
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an analytic unstable mode. Since the homogeneous collapse solution is known to suffer the 
kink instability, this result and recent numerical simulations strongly support a proposition 
that the Larson-Penston solution will be realized in astrophysical situations. It is also shown 
that the Hunter (A) solution has a single unstable mode, which implies that it is a critical 
solution associated with some critical phenomena which are analogous to those in general 
relativity. The critical exponent 7 is calculated as 7 ~ 0.10567. In contrast to the general 
relativistic case, the order parameter will be the collapsed mass. 
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APPENDIX A: ON THE PERTURBATIONS OF THE HOMOGENEOUS 

COLLAPSE SOLUTION 

Changing the position variables from r to x as 
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r = a{t)x, 



(Al) 



the time derivative at fixed r of a function / = f{t, x = r/a) is 




(A2) 



where the gradient with respect to x at fixed time is 




(A3) 



Then the equations ( |2.1|) -( ^I^ become 
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dp a dp 1 9 2 \ 

— — X ~\~ ( X pv ) 

dt a dx ax"^ dx 

dv V — ax dv dp 



0, 
GM 



dt 



dx pa dx 



a^x"^ 



0, 



dM v-dx dM 
+ 



dt 
IdM 



dx 



0, 



Aiia'^x^p. 



a dx 



(A4) 
(A5) 
(A6) 
(A7) 



We will write the perturbations of the homogeneous background as 

V = ax + u{x, t), 

p = Pb{t)[l + 6ix,t)], 
4 

M = -vrpfeaV + A(x,t), 



(A8) 
(A9) 
(AlO) 



where the expansion factor a{t) and the mass density of the homogeneous background pb{t) 
satisfy 

1 



a = aot^ pb 



(All) 



where oq is a constant. We obtain the time evolution equation for the mass density contrast 
from the equations (|A^)-( [A7| ) in the linear perturbation theory (see |^ P. 116) , 



^ + 2^— = A7TGpb6 + ^— — (x^—] 
dt"^ a dt a^ x"^ dx \ dx J 

We assume the time dependence of the perturbation as (|2.24|) , 



(A12) 



5 = t-e^W. 



(A13) 



In the homogeneous model, ([A12|) becomes 



H + -A' 

3t Vt t 



The above equations have two and only two eigenvalues a = 1 or —2/3. o 
ponds to a stable mode while cr = 1 corresoponds to a ghost mode. 



(A14) 



-2/3 correso- 
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APPENDIX B: GHOST MODE 



It is found that the system has a ghost mode, a — 1. The mode functions are given by 



Pi = P'e^ (Bl) 
Ui = U'e\ (B2) 



This mode corresponds to the following transformation: 



t^t-e, (B3) 
r^r, (B4) 



or, equivalently. 



T^T- ee", (B5) 
z ^ z + ee^z. (B6) 



This mode has no physical meaning for stability. 
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FIGURES 
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FIG. 1. U = —v/cs ioi self-similar solutions are plotted. The Hunter (A) and (C) solutions 
terminate at the second sonic point zU = —1 + z. 
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FIG. 2. zU = —zv/cs for self-similar solutions are plotted for z < 0. —z — > oo corresponds 
to the center. The Larson-Penston solution has no node which means that it is a pure collapse 
solution. Some oscillations can be seen in the Hunter (A)-(D) solutions. The number of oscillations 
is one, two, three and four for the Hunter (A), (B), (C) and (D) solution, respectively. 
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FIG. 3. P = AirGr^ p/c^ for self-similar solutions arc plotted. The Hunter (A) and (C) solutions 
terminate at the second sonic point and their density are finite there. The homogeneous collapse 
solution has a big-crunch singularity at 2; = 0. For the Larson-Penston solution, the Hunter (A) 
and (C) solutions, the density is diluted in the z > 0{t > 0) evolution. 
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FIG. 4. P = ATrGr^p/c^ for self-similar solutions are plotted for z < 0. —z^oo corresponds 
to the center. The homogeneous collapse solution has a big-crunch singularity at 2; = 0. Some 
oscillations can be seen in the Hunter (A)-(D) solutions. The number of oscillations is one, two, 
three and four for the Hunter (A), (B), (C) and (D) solution, respectively. No oscillation can be 
seen in the the Larson-Penston solution. 
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FIG. 5. z^P = AirGt^p for self-similar solutions are plotted for z < 0. The Hunter (D) solution 
has the highest central value among these solutions. 
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FIG. 6. An unstable mode funetion of the density perturbation z^5P = AirGt^dp for the Hunter 
(A) solution is plotted. A node can be seen in this mode function. It has a large amplitude near 
the center and a small amplitude with the opposite sign near the sonic point. This perturbation 
makes the concentration of the density strong or weak. 
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FIG. 7. An unstable mode function of the velocity perturbation z6U = —z5v/cs for the Hunter 
(A) solution is plotted. No node can be seen in this mode function. It is found that the "positive" 
perturbation enhances the collapse of the gas spheres, while the "negative" one promotes the gas 
to disperse away. 
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